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y^ ' We have investigated spreading of superfluid He on top of polished MgF2 and evaporated Si02 

r^ I substrates. Our results show strongly varying contact angles of - 15 mrad on the evaporated layers. 

According to our theoretical calculations, these contact angles can be explained by a spatially varying 

distribution of vortex lines, the unpinning velocity of which is inversely proportional to the liquid 

O '. depth. 

o 

(N 

Wetting phenomena in superfluid '^He have been under intense investigation during the past few years. After the 
theoretical prediction by Cheng et al. [Q, it was found out by Nacher and Dupont-Roc g] that, indeed, the alkahne 
QO ' metal Cs is not wetted by superfluid '^He. Recently, Klier et al. measured 9 = 48° for the contact angle on top of 
_1^ an evaporated Cs film using capillary rise methods. Subsequent optical experiments jjj^l gave slightly smaller values 

^ = 25 — 32°, and they also revealed strong hysteresis between advancing and receding liquid fronts. This hysteresis 

U cannot be fully explained in terms of external disorder, although roughness of the advancing contact line has been 

• found to agree with classical scaling laws [p|. 

C^ Similarly, hysteretic variation of small contact angles was observed in optical experiments where spreading on top 

G ' of a commercial antireflection coating was investigated . In this Letter we report a detailed investigation on this 

kind of contact lines and present results, both theoretical and experimental, which show that a high density of pinned 

superfluid vortices leads to apparent contact angles of a few degrees in the surface profile. In our experiments at 1 K, 

such contact angles were visible on top of evaporated Si02 layers but not on a smooth, bulk MgF2 substrate. The 

^ ■ difference follows from the larger amount of pinning sites for superfiuid vortices on the evaporated films than on a 

'~^, polished surface. Our results suggest strongly that intrinsic disorder, i.e. pinned vortices cannot be neglected when 

I ■ considering the contact line dynamics of ^He-II, e.g. on evaporated Cs films. 
^ r Vortices are always present in thin '^Hc superfiuid films; even in bulk liquid it is difficult to reach a vortex-free state 

0^ [Qjil- This is attributed to the fact that vortex core radius ag is about 1 A [g[ in ''He-II and, hence, pinning sites of 

C atomic size b can trap vortices strongly. Due to the atomic size, a large number of pinning sites can be found on the 

)^ ' surface, easily up to densities on the order of 10^^ m~^. 

Unpinning of vortices in a thickening film is caused by the interaction between neighboring vortices. A vortex 
^v ' will adjust itself to the fiow field caused by an adjacent vortex line by bending [ |lO[ . Once its radius of curvature is 
^\ , smaller than the liquid layer thickness h, the vortex becomes unpinned. Using a hemispherical pinning site of radius 

b, Schwarz calculated the unpinning velocity vq numerically llQ] . His results can be summarized by 
C^ ! 

Vq = {K4/ATrh)ln{b/ao) , (1) 

^ , where K4 = 2TTh/m4. Since the pinning potential is deep, ~ 100 K for a typical atomic site of fe = 1 nm, thermally 



iv ' activated unpinning at 1 K is effective only in the vicinity of vo ]Tl| and the geometric law of Eq. (||) is a good 
^ , approximation in our region of interest. If we assume that unpinning is caused by a single nearby vortex inducing a 
ILJ ■ velocity v = (K4/27rr), then we obtain 2h for the minimum intervortex distance within logarithmic accuracy. Hence, 
the maximum areal density no of vortices in a spatially varying film can be estimated by the relation n^ 



C^ ■ nQ0) = =r . (2) 

where K denotes the position in the plane of the substrate. 

The free energy of a superfluid film with a vortex density can be written as 
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Es= f npU{-t)dV + f adS , (3) 

Jv Js 

where Up is the particle density and a is the surface energy per unit area. The potential U{~r') denotes the energy 
per atom of mass 7714: 

C/(T^) = m4gz-^ + -^^^ln,(i?)ln(Vao) , (4) 

z-^ 47r Up 

where 7(2) is a function governing the strength of the van der Waals interaction (vdW), which decreases towards 
higher z due to the retardation effects that are known to be important for distances larger than about 5 nm ||l3|| . The 
last term denotes the kinetic energy due to vortices with superfluid density ps- The distance I is the upper cut-off 



for the vortex flow field that can also be expressed as ~ l/2dny{n). We set the logarithmic term equal to one in 
our calculation, which means that we underestimate the influence of vortices a bit. We also neglect the retardation 
effects. The vdW term in Eq. (H) is approximated by the interaction energy between two flat surfaces, v4/f 27r/i^ per 
unit area, where A is the conventional Hamaker constant. 

The minimization of Eq. (j^) depends crucially on the behavior of vortex density with increasing film thickness. 
When nv(R) is below the critical unpinning limit given by Eq. (g), the kinetic energy of vortices i?™'"*(/i) ex h. On 
the other hand, when vortices start to become unpinned and their density follows the maximum value of l/4/i^, then 
£'J°''*(/i) oc 1/h and, in fact, the gradual elimination of vortices by unpinning favors thickening of the film. This 
behavior leads to two stable minima in Es{h) as is illustrated in Fig. f without surface tension. 

The lower local minimum at film thickness hi is given by the competition between vdW forces and the elastic 
tension of vortices, which leads to the relation 



hi^ 
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(5) 



By setting n^ to its upper limit uq — 1/Ah\, Eq. (H) can be used to solve for the minimum thickness hi — ht below 
which no vortex-thinned film can exist. When h < ht, the van der Waals force wins over the vortex tension caused 
by any vortex density below the limit of Eq. (0). In other words, the vdW attraction will always limit the vortex 
density below the value of ng'"'* — l/Ahf. Using ln(Z/ao) = 1, ps = 140 kg/m^, and A = —5.9 • 10~^^J measured for 
CaFa g, we obtain ht = 12 nm and ng"* = 1.9 • 10^^ m'^. 

At a larger thickness, another minimum is obtained when the decrease of energy due to elimination of vortices is 
balanced by an increase in the gravitational energy, which yields the value 



h2 = 



K|psln(?/ao) "'^ ^ 
Wirgp 



(6) 



when rivo > 3 • 10"'^° m~^. Here we assume that there is a fluid reservoir at the level of the substrate. Note that the 
above value is much larger than the thickness of a regular vdW film because the vortex energy dies away as 1/h. A 
transition between thin (Eq. ^ and thick (Eq. |6|) films leads to an abrupt kink in the slope, reminiscent of a contact 
angle in the surface profile. Since the substrate is completely wet, we call this edge a pseudo-contact angle in order 
to make a clear distinction with a true three phase contact line. 

In our numerical calculations the one-dimensional surface profile h{x) was taken to be defined in N equidistant 
points xi, ..., xjv along the substrate. Then the minimum of the discretized free energy E(hi, ..., /iat), where hi = h(xi), 
was determined numerically with standard methods under the requirement that the total volume of the fluid remains 
constant. For the surface tension we used the zero-temperature value of cr = 375 /iJ/m^ |Tq]. 

Fig. 2 illustrates the behavior of different energy terms when a cross-over from a thin to a thick film takes place. On 
the thin film side the vortex density is taken to be n.uo = 10^'^ m^^, which corresponds to a unpinning thickness of 160 
nm. For filnrs thicker than that, ny{R) follows its maximum value uq given by Eq. (H). First, when the film thickness 
increases from 70 nm, the flow energy grows quickly (oc h), followed by a strong increase in the surface energy and 
by a rapid decrease in the vdW energy. The kinetic energy reaches its maximum at the thickness when vortices start 
to become unpinned. Above the unpinning threshold, the surface energy gradually relaxes to zero, while the vortex 
energy is reduced to a value balanced by the gravitational potential of the fluid. The inset displays the corresponding 
surface profile. An angle of 31 mrad is visible as the film breaks off from the substrate. For the pseudo-contact angle, 
we take the maximum slope of the surface with respect to the substrate. 



The calculated pseudo-contact angle 9p as a function of the vortex density Uyo is displayed in Fig. 3. Results at 
small 77.1,0 calculated on a substrate inclined by 3.5 mrad from the horizon show that the pseudo-contact angle goes 
smoothly to the value of the substrate inclination, i.e., there is no observable contact angle at small vortex densities 
UyQ ^ 10^*^ m~^. The calculations on horizontal substrate, on the other hand, indicate a discontinuous vanishing of 
Op at Uyo = 3-101° m-2. 

Our experiments were performed on a dilution refrigerator equipped with a two-beam interferometer; for details we 
refer to Ref. pq |. The investigated substrates were inclined by ^5 mrad from the horizon. We made experiments on 
three different substrates: 1) Bulk MgF2, 2) AR-coating of Hebbar-type, and 3) AR-coating of F-type 0. In our 
measurements at 1 K and below, bulk MgF2 was always covered by a smooth ^He-II film whereas substrates 2 and 3 
were not; we limit our observations to low temperatures because, in this range, the regulation of liquid level worked 
well and stable fluid fronts could be obtained. 

Wetting behavior on top of the polished MgF2 substrate is illustrated in the upper part of Fig. 4. It displays a 
smooth surface profile, deduced from the interferogram shown in the inset. The shape of the surface can be accounted 
for by the regular van der Waals film formula as can be seen from the fit to the data. This kind of a film is well in 
accordance with the measurements by Sabisky and Anderson using CaF2 and SrF2 substrates Q . The lower part of 
Fig. 4 illustrates the behavior of ^He-II on top of evaporated H ebbar-coating. A distinct pseudo-contact angle of 6 
mrad is observed in the profile. Similar surface shapes were also seen on top of the ^-coating. 

Strong hysteresis between advancing and receding front was observed on the V^-coating as well as on the Hebbar- 
coating. The contact angle was found to be zero for receding liquid whereas the advancing front displayed pseudo- 
contact angles in the range of - 15 mrad. Moreover, when shaking the cryostat, the liquid left behind tracks which 
were preferred by the next advancing front. All these features are nicely explained by a variation in the density of 
trapped vorticity in thin superfluid films. 

According to Fig. 3, the experimentally measured pseudo-contact angles of 1 - 15 mrad correspond to n„o = 
10^" — 5 ■ 10^1 ToT^. These densities are clearly larger than the values measured by Ellis and Li |]18[ who got values 
on the order of 10^ m^^ in their third sound measurements. However, all our densities are much smaller than the 
physical maximum value l.OTO^^ m~^ limited by the vdW attraction (see Eq. (||)). 

It has been suggested by Williams and Wyatt M] that surface charge might play a role in helium spreading on top 
of an AR-coating. Because of the polarizability of helium atoms, the energy density in the areas exposed to an electric 
field Eq decreases by ASnq — ^x^qEq where x is the dielectric susceptibility. If strong enough electric fields with 
sharp spatial gradients were formed, the surface profile could display distinct borders with pseudo-contact angles. 
However, the irregular behavior observed in our experiments is against the model of Williams and Wyatt: if frozen 
charge is important then a rather reproducible pattern should be seen when the liquid front is taken back and forth 
over the substrate. 

Our substrates were checked for uniformity using an atomic force microscope (AFM). The roughness was found to 
be 0.4, 1.0, and 1.0 nnirms for samples 1-3, respectively. In addition, morphology of the roughness was found to be 
different: for sample 1, the roughness was mostly caused by long, continuous scratches, while samples 2 and 3 showed 
irregular, granular-like patterns. Hence uniform, "granular" roughness on mesoscopic scale appears to favor vortex 
pinning as expected. Moreover, we used the AFM to deposit surface charge and to follow directly its decay. We found 
that the lifetime of charge on our substrates was on the order of ten minutes, which was far too short to leave any 
appreciable amount of charge behind after the mounting and cool-down procedures lasting for several hours. 

Recently, Herminghaus [g^l considered the effect of Bernoulli pressure on contact lines. His calculation basically 
applies to the spreading situation where vortex nucleation limits the superfiow and thereby governs the spreading 
dynamics. He obtained contact angles on the order of 5 degrees. However, in order to get stable liquid fronts with a 
finite contact angle in this way, a continuous heat fiux should be present in our experiments, which is not the case. 
Moreover, the pattern of the heat flux should be a rather complicated, time-dependent function of the position on the 
substrate. 

On the basis of the present experiments and calculations, we believe that the elastic energy of pinned vortices 
influences spreading of helium on top Cs [0-p| , causing extra hysteresis in the indigenous contact angle and pinning 
of the contact line. In the case of Cs, even an order of magnitude larger density of vortices is possible, since the 
evaporation of metal is done at low temperatures in order to avoid oxidation. Plenty of pinning sites were observed 
in recent flow experiments pi where solid air or solid hydrogen had been evaporated at low temperatures. 

To conclude, quantized vortices play a role in the spreading of superfluid when there are plenty of pinning sites for 
vortices. The elastic tension of vortex lines leads to thin film sections that look like non-wetted regions even though 
they are covered by *He-II. Our calculations on pseudo-contact lines at the borders of these thin film regions are 
in good agreement with experimental results obtained on evaporated antireflection coatings. Owing to unpinning of 
vortices in thick fluid layers, strong hysteresis with respect to advancing and receding liquid front is produced. This 



finding corroborates with our results, as well as recent experiments on Cs films. Hence, the dynamics and pinning of 
a superfluid contact line can be strongly infiuenced by vortices on atomically rough substrates. 
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FIG. 1. Energy per unit area vs. the thickness h of the superfluid layer. Total energy neglecting surface tension is drawn 
by the solid line. The substrate level is taken as zero for the gravitational energy, the vdW energy is zero a.t z = oo, and the 
vortex contribution is calculated for the density n„o = 10^^ m~^. For other parameter values, see text. The two minima hi 
and /i2 given by Eqs. (™ and (0), respectively, exist only when 3 • 10^" m^^ < n„o < 1-9 • 10^^ m~'^. 

FIG. 2. Energy densities and the surface profile (inset) calculated at n„o = 10^'' m~^ across the thin-to-thick film transition 
region. The solid curve is surface tension, dash-dotted is vortex energy and the dashed curve is the van der Waals interaction 
energy. The maximum slope yields the pseudo-contact angle 31 mrad. 

FIG. 3. Pseudo-contact angle 9p vs. vortex density n„o calculated for substrate inclination of (o) and 3.5 mrad (•). 

FIG. 4. Interferograms and the corresponding surface profiles measured on bulk MgF2 substrate (upper frame) and on 
-ffe66ar-coating (lower frame); the dashed lines denote the substrates inclined from the horizon. The profiles have been calculated 
along the white lines marked in the interferograms. The surface on top of MgF2 fits the standard vdW film formula illustrated 
by the solid curve. On the _ffe66ar-coating a pseudo-contact angle 9p — 6 mrad is visible. 
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